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1. Introduction 

Polymers are very long chain-like macromolecules that play an important role in a wide 
variety of physical systems. Many of the materials that we encounter in our every day lives, 
such as plastics and rubber, are essentially a mesh like structure of polymers. In biology, 
polymers also play a central role, as many biological molecules, such as DNA, have a long 
chain-like structure. In this chapter we focus on the equilibrium statistical mechanics of a 
single polymer chain immersed in a quenched random medium. For example, a chain that 
is free to move in a porous material, such as a sponge or a complicated gel network. This 
problem is relevant to a number of technologically important processes, such as filtration 
[IH 13 5 gel permeation chromatography [E2 E], and the transport of polymers through 
porous membranes [EJ. Here, we will not address these practical issues, but rather, give 
an overview of the theoretical aspects of the problem. 

In this chapter we present analytical and numerical results on the conformational statis- 
tics of polymers in quenched random media. We first consider the problem of a Gaussian 
polymer chain in a random potential characterized by short ranged correlations. This 
simplified mathematical model of a real polymer is tackled analytically using the replica 
method, and also via a mapping to the equivalent problem of a quantum particle in a 
random potential. Here, we will focus on the statistical properties of a long polymer 
chain, where the free-energy landscape is complex and possesses many meta-stable states. 
Using the path integral mapping between the partition function of a Gaussian polymer 
chain and the imaginary time Schrodinger equation, we show that the glassy phase can 
be understood in terms of the localized eigenfunctions of the Schrodinger equation with a 
random potential. Furthermore, we show that the glassy behavior can also be described 
analytically by a one-step replica-symmetry-breaking (RSB) solution. We explore the con- 
nection between the replica solution and the eigenfunctions of the Schrodinger equation, 
and show that the one-step RSB solution can be interpreted in terms of the dominance 
of localized tail states. 

We proceed to investigate the more realistic case of a chain immersed in a sea of hard 
obstacles that are randomly distributed in space. Prior to our work it was often assumed in 
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the literature that that this problem is equivalent to the random potential problem. Using 
Flory type free energy arguments we elucidate the similarities and differences between this 
case and that of a random potential with short ranged correlations. In particular, we show 
that the dependence of the polymer size on chain length can exhibit three possible scaling 
regimes depending on the system size. 

Finally, we introduce a more realistic model that includes a self-avoiding interaction 
between monomers of the polymer chain. We show that in the limit of a very long 
chain, and when the self-avoiding interaction is weak, the equilibrium chain conformation 
consists of many blobs with connecting segments. These blobs are situated in regions 
of low average potential, in the random potential case, or in a region of low density of 
obstacles in the random obstacles case. We also show that as the strength of the self- 
avoiding interaction is increased relative to the strength of the random potential, the 
polymer chain undergoes a localization-delocalization transition, where the chain is no 
longer bound to a particular region of the medium but can easily wander around under 
the influence of a small perturbation. 

2. The statistics of a Gaussian chain in a random potential 

2.1. Path integral formalism 

A polymer is a collection of molecules, called monomers, which interact with each other 
to form a long flexible chain. For example, a typical polymer like polyethylene consists 
of a chain of roughly 10 5 GHi molecules. The large number of monomers allows for a 
statistical description of macroscopic properties which are independent of details on the 
monomer scale. The simplest mathematical model of a polymer chain is referred to as the 
Gaussian chain [ 6 j . In this model the polymer is described by a position (^-dimensional 
vector R(m), where u is a continuous variable which satisfies < u < L, and which runs 
along the contour of a chain with length L. The probability of finding a conformation 
R(lt) is given by 



where M is a normalization constant, and where b is the average bond length. In this 
description the self-avoiding interaction of the chain with itself is not taken into acount, 
and in the absence of an external potential the conformation of the chain resembles the 
trajectory of a random walker with a mean step b. Thus, the average end-to-end distance 
of the chain satisfies 



The conformational statistics of a polymer chain will change if the chain is placed in a 
random environment. We can model the effect of the random environment by introducing 
an interaction energy 




(1) 



((R(L)-R(0)) 2 ) = L6 2 . 



(2) 




(3) 
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where V(R) is the potential energy of a monomer at the position R due to the environ- 
ment. The potential function V(R) will depend on the type of random medium that is 
being studied, and will be discussed in more detail later. The conformational probability 
is now 



P[R(u)} = N exp 



d 

2b 2 



—r^ / du 



dR(u) 
du 



P / du V(R(u)) 



(4) 



where we multiplied the free Gaussian conformational probability with the Boltzmann 
factor exp (—/3E int [R(u)]), with (3 = 1/kT. Given the conformational probability we can 
we write the partition sum (Green's function) for the paths of length L that go from R 
to R' as 



/•R(L)=R' 

Z(R, R'; L) — [dR(u)}exp(-(3H), 
</r(o)=r 



(5) 



where 



H 



r*T( 



M (dR(u) 



du 



+ V(R(u)) 



(6) 



and where M — dj ((3b 2 ). All the statistical properties of the polymer chain can be derived 
from this partition sum. 

In this chapter we will also consider the effects of the volume of the random medium 
i.e. the system size. In order to incorporate this effect in the partition function we include 
a harmonic term 



V(R(u))^V(R(u)) + ^(u), 



(7) 



where the coefficient fi will be a measure of the available volume to the polymer (A larger 
\i implies a smaller system volume). The confining well is also important to ensure that 
the model is well defined, since it turns out that certain equilibrium properties of the 
polymer diverge in the infinite volume limit (// — > 0). 

A quenched random medium, such as a rough surface or a frozen gel network, is a 
complex structure that can in principle be modelled by a complicated potential function 
V(R). However, we will not be interested in the physical properties of a polymer chain 
immersed in a specific environment, but rather in an ensemble of similar environments. 
Hence, we will have to specify instead the probability distribution of the random po- 
tential V(R). Here, we will consider random potentials that are taken from a Gaussian 
distribution defined by 



(V(R)} = 0, (V(R)V(R')) = / ((R - R') 2 ) . 

In particular we will consider a correlation function of the form 

9 



f ((R - R') 2 )) 



exp(-(R-R')7£ 2 ), 



(8) 



(9) 
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where g determines the strength of the disorder and the parameter £ conveniently controls 
the correlation range of the random potential. Here, we will consider only the case of 
short-range correlations, where £ is much smaller than the system size. 

All the statistical properties of the polymer will depend on the partition sum. For 
instance, the average end-to-end distance of a polymer chain that is free to move is given 
by 



/ dRdR'(R - R) 2 Z(R, R ; L) 



™ = { jitmimM) ) • (10) 

where the over-bar stands for the average of the ratio over realizations of the random 
potential. This average is referred to as a quenched average, as opposed to an annealed 
average, where the numerator and denominator are averaged independently. In some 
previous studies it has been argued that for the mean end-to-end distance, as defined in 
Eq. one can replace the quenched average by the more analytically tractable annealed 
average. However, this replacement can be justified only when the system size is strictly 
infinite, since only in that limit can the polymer sample all of space and find the most 
favorable potential well that will be similar to its environment in the annealed case. The 
main problem with this approach is that in practice we always deal with finite-size systems, 
and it is not always easy to assess how big the system size has to be so that the annealed 
average is a good approximation to the quenched average. In addition, the time it takes 
the chain to sample a large volume is exceedingly long and unreachable over a reasonable 
experimental time. 

Before ending this section we would like to point out the relationship between the 
statistical properties of a polymer chain and a variety of other physical problems. The 
partition sum, as written in Eq. (JHJ), is in the form of a path integral over all possible 
chain conformations. Since many physical problems can be formulated in terms of path 
integrals, we can map the problem of a Gaussian polymer in a random media to a wide 
variety of seemingly unrelated problems. First, we can map the partition sum of a polymer 
chain to the density matrix of a quantum particle. The mapping [[7J|H] is given by 

(3 -> 1/h, L -> ph. (11) 

The density matrix of a quantum particle at inverse temperature j3 is then related to the 
partition function of the polymer as 

p(R, R; (3) = Z(R, R; L = (3H, (3 = 1/h). (12) 

The monomer label u is now interpreted as the Trotter (imaginary) time, and M as the 
mass of the quantum particle. The density matrix is relevant to the equilibrium statistical 
mechanics of a quantum particle, such as an electron in a dirty metal. 

The polymer partition sum for a Gaussian chain in random potential can also be mapped 
into the partition sum of a flux-line in type-II superconductors in the presence of columnar 
disorder [|UJEIj- Here, R(w) is the transverse displacement of the flux-line (or vortex). 
The variable u is interpreted as the distance along the z-axis (direction of the magnetic 
field) , and the variable M corresponds to the line tension of the flux- line. Thus a flux-line 
in three dimensions propagating along the z-direction maps into a polymer in two spatial 
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Polymer chain 


Quantum particle 


Flux line 


u 

(monomer label) 


T 

(imaginary time) 


z 

(plane label) 


/? 


l/h 


(3 


L 

(chain length) 


h(5 


L z 

(distance along z-direction) 


d/(3b 2 

(6=bond length) 


m 

(mass) 


Q 

(line tension) 



Table 1 

The relationship between different physical systems. 



dimensions (its projection onto the plane), and the columnar disorder maps into point 
disorder. Although, in the case of several flux-lines there is a repulsive electromagnetic 
interaction among them which is not present for polymer chains. In Table ^ we have 
summarized the relationships between a polymer chain, a quantum particle, and a flux- 
line in a superconductor. 

Finally, we can also map the polymer problem to the problem of diffusion in a random 
catalytic environment [El^]. This process describes, for instance, auto-catalytic chem- 
ical reactions in a disordered background, or the spreading of a population with a growth 
rate that depends on local random conditions. If we make the replacements 

-j^ - D , PV(R) -> -U(R) , L^t (13) 

where D is the diffusion constant, and U(R) is the growth rate (or reaction rate) at 
position R, and t is the time. Then Z(R!, R; t) gives the concentration of the constituents 
at position R' at time t, given a delta function concentration concentrated at R at time 
t = 0. 

2.2. Flory arguments 

A remarkably fruitful approach to the statistics of polymers is via simple and intuitive 
free energy arguments. This approach was used with great success by Cates and Ball [ 
ITB*] to derive the essential scaling properties of a polymer in a random potential. Here, 
we reproduce their beautiful intuitive arguments in order to elucidate the effect of a finite 
volume on the behavior of an untethered free chain in a random potential. First, Cates 
and Ball argue that a Gaussian chain situated in an infinite random medium is always 
collapsed in the long-chain limit. Their argument goes as follows: Consider a white-noise 
random potential V(x) of zero mean whose probability distribution is 

P(V(x)) oc g' 1 ' 2 exp(-V 2 /2g). (14) 

If we now coarse-grain the medium and denote by V the average value of the potential 
over some region of volume Q, then the coarse-grained potential will have the distribution 



P n (V) oc {g/n)- 1 ' 2 e W (-nV 2 /2g). 



(15) 
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Now, consider a polymer chain situated in the random potential, and assume that it 
shrinks into a volume Q corresponding to a place where the mean potential V takes on a 
lower value than usual. In this situation the free energy of the chain is crudely estimated 
to be ( neglecting all numerical factors): 

F(Q,V) = L/R 2 + LV + nV 2 /2g. (16) 

Here, L is the length of the chain (number of monomers), R is the radius of gyration (or 
end-to-end distance) and the volume Q is related to R via Q ~ R d in <i-spatial dimensions. 
The first term on the r.h.s. is an estimate of the free energy of a long chain confined to 
a region of size R in the absence of an external potential (see e.g. [ 14 j, Eq. 1.12). The 
second term is just the potential energy of the chain in the random potential of strength V. 
The third term arises from the chance of incurring a random potential of strength V. The 
quantity In P(V) gives an associated effective entropy for the system. Minimizing this free 
energy over both V and Q determines the lowest free energy configuration. Minimizing 
with respect to V yields V = —Lg/Q. Substituting in F gives: 

This shows that for any d > 2, F — >• — oo as-R^O. Thus, the mean size of the chain is 
zero, or in the presence of a cutoff, the size of one monomer i.e. 

R ~ 1, d > 2. (18) 

For d < 2, the free energy has a minimum for 

R ~ (Lg) 1 ^- 2 ^ , d < 2 , (19) 

which in the long chain limit (L > 1/g) cuts off again at R ~ 1. These results are 
the same as those for the case of an annealed potential that is able to adjust locally to 
lower the free energy of the system. The reason is that for an infinite system containing 
a finite (even though long) chain, space can be divided into regions containing different 
realizations of the potential, and the chain can sample all of these to find an environment 
arbitrarily similar to that which would occur in the annealed situation. 

These results stand in contrast to the replica calculation of Edwards and Muthukumar 
(EM) who found that for a long chain 

Rrvg-VW, d<A (20) 

when g 2 /^~^L -> oo , whereas R ~ L 1/2 when g 2 l { ^L -> 0. Note that the result flUD 
is independent of L as opposed to Eq. (JTHJ). To reconcile the two apparently different 
results, Cates and Ball argue that the quenched case is different from the annealed case 
only for the case when the medium has a finite volume V. In a finite box, arbitrarily deep 
potential minima are not present. Instead the most negative V averaged over a region 
of volume fl V occupied by the chain, is approximately (keeping only leading terms 
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in the volume V) given by solving the equation (the l.h.s. of which represents the area 
under the tail of the distribution) 

v n 

dy P n (y) ~ - , (21) 

oo " 

which yields 



V-J'-%>. (22) 

This expression when plugged into Eq. (fT^j) leads to (note that the last term in (fTH)) just 
becomes a constant independent of R) 



When this free energy is minimized with respect to R it gives rise to 

i?~ (g\nV)- 1,{i - d \ d<4 . (24) 

Using this value the binding energy per monomer becomes 

U bmd /L = (g\nV) 2/ii - d \ (25) 

For the polymer to be localized its total binding energy has to be greater than the trans- 
lational entropy In V, which is always satisfied for 2 < d < 4 when L or V are large (and 
for d = 2 when L is large). 

2.3. The relation to the localization of a quantum particle 

In order understand the conformational statistics of a Gaussian chain in a random 
potential, we map the partition sum to an imaginary time Schrodinger equation. This 
mapping (see Ref. [|7] Eqs. (3.12)-(3.18)) is given by 



Z(R,R';L)=/ [dR(u)]exp(-pH[R(u)]) = (R\exp(-pLH)\R'), (26) 

</R(0)=R' 



where 



A =-2WW + I A2 + V{A) - (2?) 

So for a given realization of the random potential the polymer partition sum can be 
expressed as a matrix element of the imaginary-time evolution operator. The matrix 
elements can be expanded in eigenfunctions of the Hamiltonian operator to yield 

oo 

(R\exp(-/3LH)\R') = ^ exp(-/3L£ m )$ OT (R)^(R'), (28) 

m=0 
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where 

H$ m (R) = £ m $ m (R). (29) 

The Schrodinger equation with a random potential is a well known problem that has 
been intensely studied for a long time [E3 El El E] • The main property that we will use 
is that when V(R) has short range correlations (i.e. correlation length is shorter than any 
other length scale in the problem), and if the system size is infinite, then in any dimension 
all eigenstates with energy below a critical energy Em ( referred to as the mobility edge) 
are exponentially localized in the form 

$ m (R) ~exp(-|R-R m |/C). (30) 

Here R m is the localization center of the m th state, and l m is the localization length of 
that state. The localization length satisfies l/£ m = f3^j2M\E m \ for E m <C 0, i.e. deep in 
the tail region. For E > Em extended states exist when d > 2. For d = 1, 2 there is no 
mobility edge and all states are exponentially localized. The states with energies E > E M 
are called extended since they are no longer localized but are spread over a finite fraction 
of the system. Also, it is known that the eigenvalues of the localized states are discrete, 
while the eigenvalues of the extended states form a continuum. 

For finite system size, or if /x ^ in the Hamiltonian given in Eq. (|27jh the above 
discussion has to be modified. First, the eigenfunctions are always discrete in any dimen- 
sions. But even in one dimension as the energy increases the width of the localized states 
eventually becomes comparable to the system size and thus a localized particle of that 
energy can go from one end of the sample to the other. Thus the distinction between 
localized and extended states becomes blurred for a finite system at energies much above 
the ground state. Nevertheless, there will still be a qualitative difference between the low 
energy tail states and the higher energy states with large localization lengths. 

All the physical properties of the polymer chain can be expressed in terms of the 
eigenstates of Schrodinger equation. For instance we can write the end-to-end distance 
for a given realization of the random potential as 

2£ m (a m /(iRR 2 $;(R) - |J rfRR$ m (R)| 2 ) exp(-(3LE m ) 
{RUL))v = E m KI 2 exp(-/?L£ m ) ' (31) 

where a m = J c£R$. m (R), and where (-)y refers to a configurational average for the case of 
a fixed realization of random potential. When L is large enough so that (E\ — E gs )L » 1, 
where E\ is the eigenvalue of the first excited state and E gs is the ground state eigenvalue, 
then only the ground state contributes. In this case we have 

, 2 , r , x _ JrfRR 2 $ 9S (R) //dBM fl ,(R) 



(R - (L)) ^- 2 /dR<MR) ~ 2 { /dR<MR) J ' (32) 

where <3> 9S (R) is the ground state eigenfunction. It can be shown that the ground state 
wave function is positive definite and so in the large L limit (R^(L)) 1 ^ 2 can be interpreted 
as the width of the ground state eigenfunction. Assuming the ground state has the form 
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given in Eq. (f3*U)l. we can write (R F (L)) V = 2d(d+l)£^ s , where £ gs = £q is the localization 
length of the ground state. Upon averaging over all realizations of the random potential 
we get that (R F (L)) = 2d(d + l)£g S , and so the quenched average of the end-to-end 
distance, in the long chain limit, is proportional to the square of the average localization 
length of the ground state eigenfunction. 

2.4. Localized eigenstates and glassy behavior 

Using the path integral mapping we can evaluate the partition sum by solving the dis- 
cretized Schrodinger equation. In d = 1 this can be accomplished by simply diagonalizing 
a.n N x N Hamiltonian matrix, with lattice spacing A = S/N, where S is the system 
size. Details of the numerical procedure are given in [EDI- Using the lattice computation, 
we explore the connection between the eigenstates of the Schrodinger equation and the 
physical properties of the polymer chain. Here, we focus on the probability distribution 
defined as 



which can be interpreted as the probability of finding a closed polymer chain of length L 
which passes through the point R (for a given realization of the random potential). We 
consider this probability distribution since it gives the most direct connection between 
the chain properties and the eigenfunctions of the Schrodinger equation. In Fig. ^ we 
plot P(R, L) vs. R for four different chain lengths. We also include a plot of the random 
potential sample that is used. From the plot we can see clearly that as the chain length is 
increased the probability distribution tends to localize around a few valleys of the random 
potential landscape. As L is increased further there is only one peak as the chain finds 
the most favorable position. These results can be explained in terms of the eigenfunctions 
of the Schrodinger equation using the expansion 



Which shows that as L is increased the localized tail states dominate the partition sum 
until only the ground state remains when (E\ — E gs )L >> 1. 

The dominance of tail states for long polymer chains leads to large sample-to-sample 
variations in measured physical quantities. This is because the partition sum, for a given 
realization of the random potential, is dominated by a few favorable conformations which 
are strongly sample dependent. In Fig. 12.41 we plot the relative sample-to-sample fluc- 
tuations of the end-to-end distance (R F (L)), as a function of chain length. The relative 
sample-to-sample fluctuation is defined as Ap/ (R, F (L)) with 



and where the averaging is done over many realizations of the random potential. From 
the plot it is clear that as the length of the chain is increased the sample-to-sample 
fluctuations increases dramatically at a chain length L c ~ 0.5. 

Large sample-to-sample fluctuations in measured physical quantities are typical of 
glassy systems, which are a broad class of systems characterized by rugged energy land- 
scapes [Ell- This is precisely the case here, where the free energy of a long chain will 




(33) 




(34) 



in 




(35) 
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10 20 30 40 

R 



Figure 1. Plot of P(R, L) vs R for four values of L. The bottom most curve is the random 
potential sample that is used. From bottom to top we use L = .3, 1, 10, 20. The model 
parameters are M = 1/2, g = 25, (3 = 1, and ji = 0.01. We use a lattice of size S=40 
with A = 0.2. The random potential is modeled by generating N numbers {T^(«)}i=i,...,jv 
which satisfy (V € (i)^(i + 1)) oc exp(-A 2 / 2 /^ 2 ), where f = 1/^2. 



0.3 



0.2 - 



0.1 - 








1 2 3 4 5 

L 



Figure 2. Plot of Aj?/(R|,(L)) vs. L. The parameters are the same as in Fig. ^ and we 
average over 1000 random potential realizations. 
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be dominated by conformations that reside in the deep valleys of the random potential 
landscape. The main new result here is that the emergence of glassy characteristics of a 
polymer chain in random media can be traced to the dominance of localized tail states of 
the Schrodinger equation. 

2.5. Averaging over disorder: The replica trick 

The numerical results in the previous section reveal that a polymer in a random poten- 
tial can be viewed as a glassy system. In order to develop this point of view, we apply 
the replica method, which has been successfully applied in the study of glassy systems [ 
The replica approach was invented in order to compute quenched averages such as 



[R 2 F (L)), as given in Eq. (jlOj). which is difficult since it entails averaging a ratio of the 



partition sum. The starting point is the formal identity 



Hz) 



lim 



Z< 



1 



n 



(36) 



which can be used to calculate the quenched free-energy 



F = -kT ln(Z) = —kT lim 



1 



n 



(37) 



where the order of the averaging and taking the n — > limits has been interchanged, 
hopefully with no ill ramifications. The average of a term like Z n , where n is a positive 
integer, is easy to implement by creating n-identical copies of the system, referred to 
as replicas, and then averaging over the random potential. This process of averaging 
results in interactions between different replicas. The resulting analytical expression is 
then continued analytically to n = 0. 

Introducing n-copies of the system and computing Z n it is straight forward to show 
that the average end-to-end distance can be formally written as 



(R|(£)> 



lim 

n->0 



/ n: =1 n: =1 - r^z^rj, iky, v 

I nLi dR a riLi dR' a Z n ({K a }, {R'J; L) 



where 



Z n ({R a }, {R'J; L) = / X\\dRa] exp(-f3H n ) 

Jn. a {0) = Ra a=1 



(3J 



(39) 



and where 



[ L du [ L du'J2f{(Ra(n)-R b (u')) 2 ) 
1 Jo Jo nh 



(40) 



Thus, the averaged equilibrium properties of the polymer can be extracted from an n- 
body problem by taking the n — ► limit at the end. This limit has to be taken with care, 
by solving the problem analytically for general n, before taking the limit of n — *■ 0. 
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2.6. The Replica Variational Approach 

In order to compute quenched averages of the polymer chain we will have to solve 
the n— body replicated partition sum given in Eq. (|4U|) . This path integral cannot be 
evaluated analytically and a variational approach has been used in Refs. [I221IE] to make 
further progress. The procedure is to follow the work of Feynman [Ej and others 
and model H n by a solvable trial Hamiltonian h n which is determined by the stationarity 
of the variational free energy 

n(F) = (H n - h n ) hn - - In J [dR x ] • • • [dR n ] ex P (-f3h n ). (41) 



Following [ 122] we use a quadratic trial Hamiltonian of the form 

~IE f dU I rfM '5Z^( R -W- R ^K)) 2 , (42) 

^° ^° ab 

where the matrix elements p a b are the variational parameters. The physical motivation for 
this ansatz is that the replica-replica interaction in the original Hamiltonian is modeled 
by a quadratic interaction which can be different for different replica pairs. Also, the 
form of the quadratic interaction was chosen specifically [ 22] to preserve the translational 
invariance of the original Hamiltonian. Apart from the /^-dependent term, (or alternatively 
in the limit fi — > 0), the original Hamiltonian H n given in Eq. (j4*U|) is invariant under the 
transformation 

R a (u) ^R a («) + C, a = l,---,n (43) 

where C is a constant vector. This reflects the fact that in the infinite system, without 
a confining harmonic term, after averaging over the random potential the interaction is 
translational invariant. A good variational Hamiltonian must be one that can preserve this 
translational symmetry. Thus the variational ansatz must be rich enough to implement 
the relation (J43|) . The ansatz chosen by EM [E3 violated this translational invariance 
for any q ^ 0, where q denoted their single variational parameter. This imposed an 
unphysical origin on the system when none existed in the infinite volume limit that they 
considered. It was actually shown in Ref. [ l2"2~] that if the variational ansatz is rich enough 
the translation invariance actually emerges from the variational extremization even if not 
assumed explicitly for the trial Hamiltonian. 

We will now assume that in the n — > limit the matrix p a b can be parameterized 
according to the 1-step replica symmetry breaking scheme of Parisi [ 121] . In the 1-step 
RSB scheme, the matrix p a b can be parameterized as (p, p(x)) with 

p(x) = { Po 0<x<x < (44) 
^ v ' \ pi X c < X < 1 , v ' 

and where x is Parisi's replica index. See Figure® for an illustration. 

Calculating the variational free energy using Eq. (J4*T]) it could be expressed [ 122] in the 
limit n —>• as a function of the four variational parameters, i.e. F = F(ft,p ,pi,x c ). 
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Figure 3. Illustration of first step replica symmetry breaking. The limits n — > and 
m — > x c have to be taken with < x c < 1. 



The stationarity of the free energy yields four coupled non-linear equations for the four 
variational parameters. In the limit of small \i (large volume) and long chains (L large) 
these equations were solved analytically. Denoting by A the combination 

A = (jL-p+ (1 -x c )px + x c p , (45) 

we find from the analytical variational solution 

\~(g |ln M |) 4 /( 4 - d >. (46) 

From this final result we can obtain the radius of gyration which can be shown to be 
proportional to A -1 / 4 . Recalling also that | In /x| oc In V, we find 

Rg ~ A- 1 / 4 „ {g |l n/ /|)- 1/(4 - d) ~ (9 In V)- 1/(4 " d) , (47) 

a result which coincides with the prediction of Cates and Ball for the case of a finite 
volume. It could also be shown that a replica- symmetric ansatz for the matrix p yields 
a result consistent with the annealed average of the disorder where the chain collapses in 
the limit of L — > oo. 

For non-asymptotic values of \i and L the non-linear stationarity equations equations 
could be solved numerically [120] using a standard iterative method [I2H]- We found that 
for a given set of parameters there is a chain length L c (which depends on the strength 
of the disorder) such that for < L < L c there is only a replica symmetric solution. 
This is the case when the variational parameters satisfy x c = 1 and so = s\. For L > L c 
there is still a replica symmetric solution but we also find an additional replica symmetry 
breaking solution. So in this regime we find an additional solution such that < x c < 1 
and So 7^ S\. In order to decide which solution correctly describes the physics in that 
regime we compare their respective predictions to the lattice computation of (Rp(L)). 



14 



Y. Y. Goldschmidt and Y. Shiferaw 




1 2 3 4 5 

L 



Figure 4. Plot of (R F (L)) vs. L. The dotted line is generated by averaging over 10000 
samples, and error bars are found by computing the standard deviation of 10 sets of 1000 
samples. The dashed line is the RS solution, and the solid line is the RSB solution. 



In Fig. H we plot the mean squared displacement (R F (L)) vs. L for a given set of 
parameters. We plot this quantity using the lattice result, and also using the two pre- 
dictions of the variational method. Note that in the labels of the plots the average over 
the disorder is denoted by a second set of brackets rather than an over-bar. For L below 
L c w 0.73 there is only a RS solution which is very close to the lattice prediction. For 
L greater than L c the RS and RSB solutions are different and it is clear that the RSB 
solution is closer to the lattice result. We can see that the end-to-end distance saturates at 
a constant value as L increases. This behavior is correctly predicted by the RSB solution 
but not by the RS solution. 

So far we have seen that the glassy characteristics of a polymer in a short range corre- 
lated random potential is closely related to the dominance of low energy eigenfunctions. 
We also know that the variational solution possesses a RS solution for L < L c and an 
RSB solution for L > L c . It is well known that replica symmetry breaking is typically 
associated with glassy behavior, and for our model we show that the onset of RSB is 
precisely when the system begins to exhibit glassy behavior. The variational parameter 
that best reveals the transition between RS and RSB is the break point x c . If x c = 1 
then there is only an RS solution, and if < x c < 1 then that corresponds to an RSB 
solution. In Fig. |3] we plot x c vs. L using the same parameters that were used in Fig. 0] 
We can see that onset of the RSB solution is at L c w .73, after which we find that the 
break point decreases like x c oc 1/L. If we compare this result to the plot of Ap/ (H F (L)) 
vs L in Fig. 12 A\ we see that near L c « 0.5 the sample-to-sample fluctuations begin to 
rise rapidly. This result provides strong evidence that when the RSB solution is valid the 
polymer chain does indeed exhibit glassy behavior. 

2.7. Physical interpretation of the 1-step RSB solution 

In this section we study the physical interpretation of the replica symmetry breaking 
solution. Our purpose is to see if the underlying physical picture predicted by 1-step 
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Figure 5. Plot of x c vs. L. The parameters are the same as those used in Fig. 0J 



RSB is indeed consistent with the presence of exponentially localized eigenstates. The 
following analysis is valid for a very long polymer (large L) when the system becomes 
glassy. We begin by evaluating the replicated partition sum defined as 



Z n ({R a }) = / Y[[dR a ]exp(-(3h n ), 

-/R a (0) = Ra a=1 



(48) 



where h n is the quadratic trial Hamiltonian in Eq. (|42j) . Since h n is quadratic the path 
integrals can be evaluated analytically and the final result can we written in the form 



Z n ({R a }) = const, x exp (-^^QjR-a • 



(49) 



The details of this calculation along with the relationship between the matrices Q a b and 
p a b are given [EDI- Now, since p a b was parameterized according to the 1-step RSB scheme, 
it implies that Q a b can also be parameterized in the same way. 

Mezard and Parisi [ El] discuss the interpretation of a representation of the form (}4T?|) 
for the case of directed polymers. In particular they show how to deduce the structure of 
the probability distribution 



Py(R) = Z V (R,L)/ J dRZ v (R,L), 



(50) 



which is the probability of finding a polymer loop that passes through R for a given 
realization (which we denote by V) of the the random potential. Here Zy{R, L) is just 
the partition sum Z(R, R;L) as given in Eq. (J3j). This probability is related to the 
replicated partition function given in Eq. by 



(51) 



P V (R) = lim / dR 2 ■■■dR n Z n ({R a }) Rl=R . 

n— +0 
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Mezard and Parisi's analysis has to be adapted for the case of real (non-directed) polymers 
of length L in a random potential which is independent of time. The changes will be 
pointed out below. 

If Q a b is parameterized by {q, q(x)} such that 

q(x) = { % X<Xc (52) 

one proceeds to obtain Py(R) by the following procedure: 

1. For each sample ( a realization of the random potential) generate a random variable 
R which is picked from the distribution 

P(R ) = - * exp . (53) 



2. Consider a set of "states" labeled by the index a whose physical meaning will be 
elucidated shortly. Each of these states is characterized by a weight W a and a position 
variable R Q . Given Rq, the variables R a are an infinite set of uncorrelated random 
variables distributed according to 

V(K U R 2 , ) = J] 1 exp f- ( ^~ R °n . (54) 

The distribution of weights will be discussed below. 

3. Given these "states" for a given sample, The probability distribution Py(R) for that 
sample has the form 

VV } ^ Q (27r(g- gi )) d/2 P V 2(g- gi ) J 1 } 

The weights W a are given in terms of some "free energy" variables f a : 

W a = ^ ■ (56) 

These free energy variables are chosen from an exponential distribution 

P[f a ]^exp(x c Pf a ) 9(f-f), (57) 

where / is an upper cutoff. 

What is the meaning of these variables in the present case? To determine the weights 
W a we compare Eq. (|55j) to the eigenfunction expansion given in Eq. (|34|). From Eq. 
together with Eq. (|5Uj) it becomes clear that 

P y (R) = ^A a |$ a (R)| 2 , (58) 

a 

where 

exp(-/3LE a ) , . 

> exp -/3LE 7 
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Comparing Eq. (J55j) and Eq. (foTlj) it becomes obvious that W a = A a and 

^(R)ocexp f- S^j^ V (60) 



2(q-qi) 

Hence, the "states" labeled by a are in our case the actual eigenstates of the imaginary 
time Schrodinger equation. These are localized tail states centered at position R a with 
an associated "weight" W a . Thus the 1-step RSB solution approximates the tail states 
by a fixed Gaussian form. 

The width of these Gaussians, denoted by wo, give an estimate of the size of the polymer 
chain i.e. wq ~ Rf- In the limit of large L and small /i it can be shown to leading order 
[E2] that w 2 ~ d/(2{3^/\M), where 

A = (^^T (^A/) <1+d,/<< -* |h^l)^<- rf > . (61) 

So in terms of the disorder strength and the system size the chain size scales like 

Rf oc (g\ ln/i|)- 1/(4 - d) oc (g In y)- 1/(4 " d) . (62) 

It should be emphasized that the subtle dependence on the volume of the system is a 
direct consequence of replica symmetry breaking. In fact, as shown in Fig. Elthe replica 
symmetric solution does not correctly describe the size of the polymer chain, since it fails 
to capture the dominance of localized tail states. 

We now consider the distribution P(R a ) given in Eq. (j54j). This is just the distribution 
for the localization centers R Q for a given value of Ro- Hence, we can calculate the average 
distance between the localized states for a given sample. We find that the width w of the 
Gaussian P(R a ) satisfies w 2 = d(q\ — q ). For small \i and large L it is straight forward 
to show that that w 2 ~ d/((3u,Lx c ), where the break point x c is given by 

\ -V(4-cQ 

L V(2^p' 

2.8. Density of states and the 1-step RSB solution 

To develop the analogy further we first notice that the free energies f a are equal to LE Q . 
This make sense if we think of \E a \ as representing the binding energy per monomer, 
and thus f a = LE a represent the total energy of the chain. These arguments lead us 
to expect that within the 1-step RSB scheme, the energy variables E a are independent 
random variables taken from an exponential distribution: 

P[E a ]^e^ Lx ^e(E-E a ) , (64) 

with E being some energy scale determined by the upper cutoff of the tail region. We will 
now argue that the distribution given above is just the expected distribution of ground- 
state energies i.e. the probability of finding the lowest energy level to have energy E. 
We first review some very basic results of extreme value statistics as presented in Ref. [ 
Given K independent and identically distributed random variables E iy pulled from a 
distribution of the form 

P(£) = A eX p(-5|£f), (65) 



x c = \{ f—g 2 f3 d+i M d |ln/i| d - 2 ) . (63) 
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the probability that the lowest of the K energies is E (for E — > — oo and K — > oo) is 
given by 

P(E) oc exp [B5\E c \ s - l E] , (66) 

where 

* = -(T)"". (67) 

The value of E c , the lowest energy expected to be attained in K trials, is easily obtained 
from 

E c 

dEP(E)~l/K. (68) 

-oo 

The reason why we chose a distribution of the form given in Eq. (JfiHj) is that in d = 1 
the probability P(E) is known to have that form exactly for the case of delta correlated 
random potentials (see [I2ZI)- For d > 1 Lifshits [ ^1 argued that the form given by 
Eq. is also valid. Our goal now is to see if the distribution Eq. (JMj) derived using the 
1-step RSB solution is indeed consistent with the distribution Eq. predicted using 
extreme value statistics. 

Comparing Eq. (J64)) and Eq. (J66)) we find that for consistency the break point should 
satisfy 

x c = ±B 1 ' s Qog{K))t'- 1 V s . (69) 

Notice that the 1/L behavior of x c is exactly the same as was found analytically for large 
L in Ref. [122] and numerically for any L > L c in the present work. We can go further 
by using the fact that the number of energy levels K, within a fixed energy interval is 
directly proportional to the system size, which in our formulation is effectively determined 
by fi. Assuming log(if) oc | log(/x)| and comparing to the approximate solution for x c in 
Eq. we find that 5 = (4 — d)/2 and B oc 1/g. Now P(E) is just proportional to the 
density of states p(E), and it is known exactly in one dimension. Indeed, when d = 1, 
5 = 3/2 and B oc 1/g. For 2 < d < 4, 5 agrees with the result derived by Lifshits [IT9"]. 
Hence, the exponent 5 and the disorder dependence of B is correctly predicted by the 
1-step RSB solution. 

The above results show that the 1-step RSB solution correctly predicts some important 
features of the eigenvalue distribution. More importantly, we have shown that the 1-step 
RSB solution can be interpreted in terms of the eigenstates of the Schrodinger equation 
with a random potential. However, there are differences and these reveal the limitations 
of the 1-step RSB solution. For example, all the localized states are approximated by the 
same Gaussian profile when in fact the localization lengths should increase with energy. 

3. Localization of polymers in a medium with fixed random obstacles 



In this section we discuss the static properties of a Gaussian polymer chain without 
excluded volume interactions that is confined to a medium populated with quenched 
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random obstacles. It is important to distinguish between the following two important 
cases that have been discussed in the literature: 

1. A Gaussian random potential with short range correlations. 

2. Random obstacles which prevent the chain from visiting certain sites. 
Numerical simulations performed in three dimensions were restricted, to our knowledge, 

only to the case of random obstacles [123123 ED]- On the other hand, extensive analytical 
work using the replica variational approach and Flory type free energy arguments, has 
been done for the case of a Gaussian random potential [ H3J EH I2D1 1221 HI] ■ The case of a 
bounded (saturated) random potential was also addressed in Ref. [H3]. It was not clear to 
us to what extent these theoretical investigations could be applied to the case of infinitely 
strong random obstacles placed randomly in the medium, as simulated numerically. This 
motivated us to investigate this problem in detail [EI]- The results indicated that only 
in a special small value of the embedding volume) the two problems mentioned 

above are similar, but otherwise they are quite different. 

We will assume that the obstacles are infinitely strong-they totally exclude the chain 
from visiting a given site occupied by an obstacle. Each obstacle is taken to be a block of 
volume a d , where d is the number of spatial dimensions and where a is the linear dimension 
of the block. We take for simplicity the polymer bond length b to be approximately equal 
to a. Thus, a will be the small length scale in the problem, and we will measure all 
distances in units of a. The obstacles are placed on the sites of a cubic lattice with lattice 
spacing a. We denote by x the probability that any given lattice site is occupied by an 
obstacle (block). Our main results will concern the case of small x, in particular x < x c , 
where x c refers to the percolation threshold (x c = 0.3116 for a cubic lattice in d = 3), but 
we will also comment on the case of a larger concentration of obstacles. We denote by V 
the total volume of the system. 

Assume that the chain is occupying a spherical region (lacuna) of volume Q ~ R d . In 
this region the actual volume fraction of obstacles will be denoted by x. It is crucial to 
realize that although the average number of obstacles per site is fixed by x, the actual 
number of obstacle in a small region of volume ~ R d is a fluctuating quantity which occurs 
with probability b(R d x; R d ,x), where 



denotes the binomial probability distribution. 

In the limit where the system has infinite volume V the free energy for the chain is 
given by 



All these terms originate from entropy F = —TS where for simplicity we take T = 1 
since the temperature does not play a significant role here with respect to the results. 
The first term originates from the entropy of a free chain in <i-dimensions where z is the 
coordination number which for a cubic lattice is equal to 2d. The second term originates 
from the entropy of confinement in a cavity of radius R. The third term is the entropy 
loss due to obstacles. This linear dependence is justified in Ref. [EI] and in the Appendix 




(70) 




(71) 
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therein. The forth term represents an entropy given by the logarithm of the probability to 
have a region of size Q with fix obstacles. This free energy, valid for V — > oo, is called the 
annealed free energy since when the polymer can sample the entire space it is the same as 
the random potential adjusting itself to the polymer configuration. The free energy has 
to be minimized (the entropy maximized) with respect to R and x. The most favorable 
value of x is 0. Since 6(0; R d , x) — (1 — x) R , we find 

F{R) = -L hx{z) + ^-R d ln(l - x). (72) 

This free energy has now to be minimized with respect to R to yield 
/ I \ V(«*+2) 

Rm.annealed ~ I 7~i 7~t \T I ■ (' ^) 

\\\n(\-x)\J 

Thus the size of the chain grows with L, but with an exponent smaller than 1/2, the free 
chain exponent. 

So far we discussed the case of an infinite volume V. In a finite volume we find that the 
so called quenched and annealed case differ, at least when the volume is not too big. We 
actually find that there are three regions as a function of the size of the system volume V. 
First, if V < Vi ~ exp(x~( d ~ 2 '' 2 / '(1 — x)), it is unlikely for a chain of volume fl ~ R d to 
find a region which is totally free of obstacles. Thus x does not vanish in this regime. To 
proceed further we must use an approximation to the binomial distribution b(flx;fl,x). 

If fl is not too small we can approximate the binomial distribution by a normal distri- 
bution 



b(flx; fl, x) « (2^(1 - x)Y 1/2 exp f-^L-E^ J . (74) 

V 2x(l — x) J 

This approximation is good provided fix ^> 1 and Q(l — x) ^> 1. We verified that these 
conditions are indeed met in our case when x is small. 

In a finite volume V, the lowest expected value of x, to be denoted by x m , can be found 
from the tail of the distribution 

( V{x-x) 2 \ Q 
dxeXP { 2^J"V' (75) 



which gives 



xy In V 



R d ' 

where we put y = 1 — x. The free energy becomes 



(76) 



Fj(R) = -L \n(z) + ± + Lx-L^I (77) 

The last term in the annealed free-energy is missing since it is negligible for large L when 
R is independent of L. Minimizing F(R) with respect to R we find 

R mI ~ (xylnVy 1/i4 - d) (78) 
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and 

x mI = x - (xylnV) 2/(4 - d) . (79) 

The result for the radius of gyration of the chain, as represented by R m i is the same 
result as for the case of the Gaussian distributed random potential, but with the strength 
g replaced by x(l — x). The polymer in this case is localized and its size is independent 
of L for large L. 

As V grows R m decreases until eventually x m vanishes. This happens when V = Vi — 
exp(x^ ( - d ^ 2 ^ 2 y^ 1 ). For V > Vi, R m is no longer given by R m i, but rather by the solution 
of x m ii = 0. It is the largest region free of obstacles expected to be found in a volume 
V. Rather than using the normal approximation we can estimate R m directly from the 
relation 

(l-x) n ~Q/V, (80) 
with fl ~ R^. Solving for R m we obtain 
( InV \ 1/d 

RmII ~{\W^T\) • (81) 

The polymer is still localized but the dependence on x and on In V has changed. In this 
region which we call region II the free energy is given by 

(l v \ ~ 2 ^ 

where some undetermined constant a is introduced for later convenience. 

As V grows in region II, Rmii continues to grow until it reaches the annealed value 
given above. This happens when 

V = V 2 ~ exp ( a ?/{M) L d/(d+2)) (83) 

to leading order in x, which is enormous for large L. For V > V2 we have the third region 
in which R mIII = R m ,anneaied and it grows like L l / {d+2 \ 

We can thus summarize the behavior of the end-to-end distance as the function of the 
system's volume as follows: 

Region I V < Vi ~ exp(2r (d ~ 2)/2 ) 



R mI ~ (xlnV)" 1/(4 - d) (84) 

Region II Vi < V < V 2 exp (^2/(^+2)^/^+2)) 

InV ^ 1/d 



^^NT? 1 • (85) 

region III V 2 < V 
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The behavior in region II can be deduced from known results of the density of states 
for a quantum particle in the presence of obstacles (repulsive impurities). In that case [ 
HHj the density of states is given by (when the obstacles are placed on a lattice) 

p(E) ~ exp(-c| ln(l - x)\E- d/2 ), E > (87) 

with c being some dimension dependent constant and x is the density of impurities. Note 
that p(E) vanishes for E < 0. We can estimate the lowest energy in a finite volume V 
from the integral 

E c 

dEp(E) ~ 1/V, (88) 



o 



and find 

„ / InV \~ 2/d , . 

ft ~(iMT^i) ' (89) 

and thus the localization length is given by 

£ „ \E I-V2 ^ f 1^ ^ 7 (go) 

We now make some remarks on the validity of the spherical droplet approximation. The 
shape of a long polymer chain is determined by the regions of the random medium that 
have a lower than average number of obstacles. For V > Vi these regions are essentially 
free of obstacles. The probability of finding such empty regions depends only on its volume 
and not its shape. However given regions of varying shapes and equal volumes, it will be 
entropically more favorable for a long polymer chain to reside in a region whose shape is 
closest to a sphere. This is because the confinement entropy is maximized for a sphere 
over other shapes of the same volume. The argument is equivalent to that proposed by 
Lifshits [ in the context of electron localization and is shown rigorously by Luttinger 
[ I33j . For V < Vi the relevant regions contain a small number of obstacles but we believe 
that the same argument should roughly hold and deviations from a spherical shape will 
be small or irrelevant. 

We have compared our analytical results with numerical simulations performed by 
Dayantis et al. [EDj, and also comment on the relation to earlier simulations done by 
Baumgartner and Muthukumar [I2H1- Dayantis et al. carried out simulations of free 
chains (random-flight walks) confined to cubes of various linear dimensions 6 — 20, in 
units of the lattice constant. These chains can intersect freely and lie on a cubic lattice. 
They introduced random obstacles with concentrations x = 0, 0.1, 0.2 and 0.3. The length 
of the chains vary between 18 — 98 steps. They also simulated self- avoiding chains that 
we will not discuss here. They measured the quenched entropy, the end-to-end distance, 
and also the radius of gyration which is a closely related quantity. Unfortunately, these 
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Figure 6. A plot of ln(—S/L + In 6) vs. ln(lnV/| ln(l — x)\). The labels are marked 
according to the chain length. 



authors did not have a theoretical framework to analyze their data, and thus could not 
make it collapse in any meaningful way. We show below how it is possible to fit the data 
nicely to our analytical results. 

Even for x — 0.1, the value that we get for Vi is about 33 which is an order of magnitude 
smaller than the the smallest volume used in their simulation, which is 216 for a cube of 
side 6. Hence we expect to be in region II. To check the agreement with our analytical 
results we show in Figure |H1 a plot of ln(— S/L + In 6) vs. ln(lnV/| ln(l — x)\) where S 
is the entropy measured in the simulations and V = B 3 for a box of side B. Recall that 
F = —S and Eq. (JH2~|) predicts a straight line with slope —2/3. The best fit is obtained 
for a slope of —0.72 ±0.05, which is in excellent agreement with our analytical results in 
region II. 

In order to analyze the simulation results for the end-to-end distance and radius of 
gyration we have to introduce some additional compensation for the results obtained in 
the previous section. First we must realize that Eq. (JHTj) is valid only when the number 
of steps (monomers) is very large. In the simulations they used chains of varying lengths 
whose size did not yet reach asymptotia. Hence, we introduce a correction factor 

R m (L) = R m {l - exp(-L/iO V2 = RmfiiVL/Rm) , (91) 

which interpolates between the size of a free chain as L — > and the value of R m from 
Eq. (jHU) as L -> oo. 

The second correction we have to implement arises when the expected value of the chain 
is not much smaller than the size of the confining box. Even for a free chain confined to 



24 



Y. Y. Goldschmidt and Y. Shiferaw 



a box of side B with no obstacles present, the end-to-end distance is not simply R = L 1 / 2 
for L 1 / 2 < B and R = B for larger L. We have to take into account the fact that the 
length of the chain has a Gaussian distribution about its expected value, and the tail of 
the Gaussian is cut off by the presence of the box (this is for the absorbing boundary 
conditions that is used in the simulations). Thus, for the case of no obstacles (x — 0), 
The measured end-to-end distance should approximately be 

R 2 C = J B JRR 2 eM-^)/ J* dR exp(-^), (92) 

which gives R c = y/Lf 2 {B/y/L) with 

^^{'-^^m^-^f 2 - (93) 

This indeed gives good agreement with the measured values in the no obstacle case. For 
the obstacle case we thus have to introduce these two corrections in subsequent order: 

Rm,corrected = Rmfl{^L/Rm)f2 ( ~ r— ) , (94) 

\R m ji[yL/R m )/ 

where R m = R m ii as given by Eq. (JHTj) . (A constant of proportionality of 1.8 has been 
introduced on the rhs of Eq. (JHTj) to obtain a good fit). In Figure 13 we show a comparison 
of the simulation results for the end-to-end distance with the calculated results as given 
by Eq. (J81j) and Eq. (|94j) . The agreement seems remarkable, since all the data collapses 
to a straight line with a slope close to 1. 

Dayantis et al. emphasize that they did not consider concentrations of obstacles above 
the percolation threshold, which is at x c = 0.3116 for a simple cubic lattice. The reason 
is that above the percolation threshold the medium of random obstacles begins to form 
disconnected islands free of obstacles. Thus, in their simulation the polymer chain will 
only sample a limited fraction of the volume available. What happens is that effectively 
the volume available for the chain is not the total volume of the cube but rather the volume 
of the disconnected region it occupies. For most realizations of the random medium this 
effective volume will be smaller than the value Vi, which is the limit of region I of the 
last section. In that case one expects the end-to-end distance to scale like x _1 as given 
in Eq. (J78)) instead of like x -1 / 3 as given by Eq. (J81j) . Baumgartner and Muthukumar's 
simulation was for both below the percolation threshold and also above it (x = 0.4 and 
0.5). However, they only estimate the exponent above the percolation threshold, and 
find it to be about —1. They do not estimate the exponent for x below the percolation 
threshold, which appears from their data to scale with a much smaller exponent. Thus, 
it seems likely that the reason these authors report a behavior corresponding to region I, 
even though their box is quite large, is because the effective volume is small for the cases 
for which they exceed the percolation threshold. 
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Figure 7. A plot of the observed vs. calculated end-to-end distance 



4. Localization and derealization of polymers with self-avoiding interaction 
in the presence of disorder 

In the previous sections an ideal (Gaussian) chain has been used, which corresponds 
approximately to the experimental situation at the so called G-temperature when the 
solvent effectively screens the self-avoiding interaction of the chain. In some early papers 
[ 123 EH EH] there was an attempt to include the effects of the self-avoiding interaction 
of the polymer. These attempts were far from complete. For example in Ref. [ it 
was assumed that the conformation of the polymer consists of one spherical blob and it 
was argued that a quenched random potential is irrelevant for a very long chain when 
a self-avoiding interaction is present. In Ref. [EH] analytical results were obtained for 
annealed disorder, and simulations were performed for strictly self-avoiding walks. Ref. [ 
|2*%] presents numerical evidence for a size transition of the polymer as a function of the 
relative strength of the disorder and the self-avoiding interaction. The simulations were 
carried out for a random distribution of hard obstacles with a concentration exceeding the 
percolation threshold. In a recent paper [ES] we tried to shed more light on this important 
problem. We have made use mainly of Flory-type arguments, and considered both the 
case of a Gaussian random potential and the case of randomly placed obstacles. Note 
that a polymer with self-avoiding interactions cannot be mapped into a quantum particle 
at a finite temperature in a simple manner, because for a quantum particle there is no 
impediment to return at a later time (or Trotter time) to a position it visited previously. 

An important point to keep in mind is the strength of the excluded volume interaction. 
If one considers a strictly self-avoiding walk on a lattice (SAW) corresponding to a non- 
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self-intersecting chain, then the strength of the Edwards parameter v [ |H| is fixed at 
0(1) x a 3 where a is the step size (or monomer size) and depends only on the type of 
lattice. On the other hand one can consider a Domb- Joyce model [ where there is a 
finite penalty for self overlapping of polymer segments, and then the strength of v can be 
varied substantially and reduced continuously to zero. The interplay between the strength 
of the self-avoiding interaction and the strength of the disorder can then be investigated 
to a larger extent. Experimentally the Edwards parameter is given approximately [E] by 
v = a 3 (1 — 2%), where \ is the Flory interaction parameter, which depends on the chemical 
properties of the polymer and the solvent, and on the temperature (and pressure). It takes 
the value 1/2 at the G-point. The case x — corresponds to a solvent that is very similar 
to the monomer. In general good solvents have low \ whereas poor solvents have high 
X resulting in v being negative. In the following we will restrict ourselves to the case of 
positive v, which leads to the more interesting and non-trivial results. 

We now revisit the meaning of the term localization as applied to polymers in a random 
medium. Although some authors connect the compact size of the chain when L — > oo 
with the notion of localization, this is actually not so. The compact size should be viewed 
as a separate feature from the notion of localization. Recall that for a Gaussian chain in 
an uncorrelated Gaussian random potential of variance g the chain has typical size 

R F oc (3lnV)- 1/(4 - d) , (95) 

and the binding energy per monomer is given approximately by 

U bind /L~-(g\nV) 2 ^- d l (96) 

To insure localization, the binding energy of a chain Uund has to exceed the translational 
entropy In V. From Eq. (|96j) this amounts to the condition 

InV < L(g\nV) 2/{4 - d \ (97) 

which holds for any 2 < d < 4 when L is large enough (for 2 < d < 4 and any fixed L, the 
condition can be satisfied for large enough V) . This condition assures that the polymer 
will stay confined at a given location and will not, under a some small perturbation move 
to a different location. Thus repeating an experiment or a simulation with the same fixed 
realization of the disorder, but with different initial conditions, will result in finding the 
polymer situated at the same region of the sample as in a previous experiment, provided of 
course one waits enough time (which can be enormous) for the system to reach equilibrium. 
We observe that this condition is satisfied for large enough L provided the binding energy 
per monomer is positive. Another interpretation of the inequality given above in the 
context of equilibrium statistical mechanics is that the partition sum is dominated by 
the term involving the ground state as opposed to the contribution of the multitude of 
positive energy extended states. The contribution of these states is proportional to the 
volume of the system and thus the inequality above results from the condition 

ex V (-LE ) > V. (98) 

What we will see in the following sections is that in the presence of a self-avoiding in- 
teraction, a localization-delocalization transition occurs when varying the strength of the 
the self-avoiding interaction for a fixed amount of disorder or alternatively upon varying 
the strength of the disorder for a fixed value of the self-avoiding interaction. 
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4.1. A self avoiding chain in a random potential 

Consider first the case of a random potential with a Gaussian distribution. For sim- 
plicity, the discussion in the rest of this section be limited to three spatial dimensions 
(d — 3). Recall that in the case when there is no self-avoiding interactions the optimal 
size of a chain R m is found by minimizing the free energy F in Eq. 1221 This yields 

Rm ~ i ~~7T7T = T^ly (") 
g\og{V) G(V) 

where we defined the volume dependent disorder strength by G(V) = gln(V). Substituting 
this result in F we obtain 

-G(V) 2 L. (100) 



in 



We see that —F m /L is the binding energy per monomer, and it is strictly positive, so the 
polymer is localized. In what follows we will assume that g is small enough so that G ^ 1 
for the given system volume, hence R m ^> 1, and the chain is not totally collapsed unless 
V -> oo. 

We now add a self avoiding interaction and assume first that it is small, i.e. v << g, or 
at least v < g. If the chain is still localized in the same well, which we will see momentarily 
not to hold when L is large, then 




flflog(V) , vL 



2 



+ — ^-G(VyL + vG(VyL\ (101) 

Here, besides assuming that v is small we assume for the moment that L is not too big 
so the last term in the free energy, resulting from the self-avoiding interaction, is small 
enough so one does not have to take into account the change in R m due to the presence 
of v. If we plot F vs. L, we see that it is lowest when 

L = L m = -— . (102) 



Thus if vg <C 1, we have L m ^> 1, and 
F -_?(V) 



(103) 



For L = 2L m the free energy vanishes and for larger L it eventually increases fast like 
L 2 . We can now verify that if L does not exceed L m then the approximation used above, 
assuming that i? m does not change appreciably from its ideal chain value, is justified. If 
we differentiate the above F in Eq. (jlOlj) with respect to R m we find 

again omitting constants of order unity. Thus the correction vGL evaluated at L — L m 
is of order unity, and we can still use the value R m ~ 1/G. 
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Figure 8. Conformation of a chain consisting of many "blobs" with connecting segments 
as L increases. Initially the first blob "overflows" and there is "hopping" to another deep 
minimum. There is a competition between the length of the connecting segment and the 
probability to find a deep minimum nearby. 



For larger L the approximation seems to break down, but fortunately what happens is 
that since the free energy increases when L exceeds L m , it is energetically favorable for 
part of the chain to jump into a distant well. Even though there is a cost for the polymer 
segment between the wells one still gains in the overall free energy from the binding energy 
in the wells. Thus the picture that emerges is that as L increases, the chain divides itself 
into separate blobs with connecting segments. In each blob the number of monomer does 
not exceed L m , which is the optimal value for that well. The idea is depicted in Figure |H1 

To be more specific we will now construct a model for the free energy of the chain. 
The first blob will be located in the deepest minimum in the total volume, whose depth 
is roughly given by — G(V) 2 per monomer, where V is the total volume of the system. 
Subsequent blobs will reside in the most favorable well within a range Y, which has to 
be taken self-consistently as the length of the jump. Thus within a range Y the chain is 
likely to find a potential minimum of depth —G(Y 3 ) 2 . The farther you jump, the deeper 
well you are likely to find. Thus assuming for simplicity that the jumps are roughly of 
equal size, and there are K blobs in addition to the initial blob, the free energy of the 
chain will be given roughly by 



F(w,m,Y;L) 



G(V) 




vm 



- wG{Y 3 ) 2 + vw 2 G(Y 3 ) 3 , (105) 
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with 

K = L = — i—, G(Y 3 ) = 3gHY ). (106) 

w + m 2vG(V) 

We defined w to be the number of monomers in each blob, and m to be the number of 
monomers in each connecting segment. The term Y 2 /m results from the "stretching" 
entropy of the segment and m/Y 2 from the confinement entropy. The term vm 2 /Y 3 
represents the self-avoiding interaction for the connecting segments. Lq is the number of 
monomers in the initial blob whose free energy was taken care of separately. It is evident 
that when L is very large we can neglect the free energy of the first blob and also take 
K m L/(w + m). Thus we find for the free energy per monomer 

,-i x f\ F(w,m,Y;L) 1 (Y 2 m vm 2 2^/^s\s 
f(w, m, Y) = v — — « h — + — =- - wG(Y 3 ) 2 + vw 2 G(Y 3 ) 3 



w + m \m Y 2 Y 3 

(107) 

This function has to be minimized with respect to w, m and Y to find the parameters 
giving rise to its lowest value. For the connecting pieces of the chain we did not include 
a contribution from the random potential since it is expected to average out to zero for 
these parts. To gain some feeling into the behavior of this function and the values of the 
parameters which minimize it, we display in Table I the value of the parameters and free 
energy per monomer for g = 0.05 and various values of v, as obtained from a minimization 
procedure. The derealization transition is the point where / changes sign from negative 
to positive, as discussed earlier. Actually, to be more precise, the derealization transition 
occurs when / = — (lnV)/L for finite L, when the translational entropy starts to exceed 
the binding energy. In the limit of large L we can say that the transition is at / = 0. 
We observe that the derealization transition occurs at v — 0.0478 which is close to the 
value of g. We also observe that m <C w for v g and m ^> w near the transition. 
Also for small v, m ~ Y, whereas near the transition m ~ Y 2 . If we compare the value 
of w from Table I with the value of L m = l/(2vG(Y 3 )) we find that w is smaller than 
L m in the entire range. The ratio w/L m varies from ~ 0.4 to 1 as v changes from 10 -5 to 
0.048. Thus the assumption we have made previously concerning this ratio is justified a 
posteriori. 

Luckily it was possible to solve the minimization equations analytically almost entirely 
in both the limits v g, and near the transition when / « 0. Details of the solutions 
are given in the Appendix of Ref. [155]. Here we only display the results: 

A. The case v <C g. 

The parameters are given by 

Y = -(lnF-n-^Hny + 3)" 3/2 , (108) 
v 

m= 3^T7 (lnl '" irl(lnl ' + 3r1 ' (109) 

W= 3^Y (,nY + 3rl (U0) 
/ = -9/(lnF) 2 (lnr - l)(lnY + 3)" 1 . (Ill) 
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V 


Y 


m 


w 


f 


0.00001 


2206 


2413 


16178 


-0.835249 


0.0001 


346 


534 


2580 


-0.421023 


0.001 


60 


148 


461 


-0.164673 


0.01 


13.7 


66 


97 


-0.0370883 


0.02 


11.1 


69.2 


60.8 


-0.0156444 


0.03 


12.8 


105.5 


41.9 


-0.0056033 


0.04 


22 


302.6 


26.8 


-0.0009992 


0.045 


35.6 


720.1 


20.7 


-0.0001776 


0.0478 


69.4 


2342.4 


16.5 





0.048 


77 


2800 


16 


0.00000725 



Table 2 

Parameters and free energy for the case g=0.05 



The first equation can be easily solved numerically for Y for a given value of v and the 
result substituted in the other equations. Very good agreement is achieved with Table I 
for small values of v. 

B. Solution near the derealization transition. 

Let us define the parameter 

Av 

«=— . 112 
In terms of this parameter we have 

Y = L(i±m?i (113) 



W 8v 2 \nY 2vG{Y 3 )' 
The transition point is obtained by solving the equation 



! + « ■ '(i-2h.r (1+ £^||=a me) 



« l + vi + k 2 k \ \ V3gK 

Once the solution k c is determined for a given g, then the transition point v c is determined 
from v c = 3g/€ c /4. For g in the range 0.01-0.2 we find that k c is a number of order unity 
(varies from 1.7 to 0.96 as g changes in that range). This means that v c is quite close to 
g. Once v c is known, all the parameters Y,m and w at the transition are determined by 
the solution above. For g = 0.05 we get k c = 1.2744, and thus v c = 0.0478 in excellent 
agreement with the minimization results from Table El 
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Figure 9. A typical chain conformation when v ^> g. The dark regions are regions of low 
average potential. Only short segments of the chain are situated in these regions. 



For v > v c the chain is delocalized. The above expression for the free energy may no 
longer be accurate, but the general picture is clear. There will be very few monomers in 
the low regions of the potential, and the chain will behave very much like an ordinary 
chain with a self-avoiding interaction in the absence of a random potential. Any little 
perturbation can cause the chain to move to a different location in the medium (see Fig. 

©)■ 

For small values of v when the chain is localized, we still expect its size to grow like 
that of a self-avoiding walk. Thus we expect roughly 

R g ~ YK 0S , (117) 

since Y is the step size, and K is the number of steps. But since K « L/(jn + w) for 
large L, we find 

R g ~ -. — -—-tL®'®. (118) 

Thus the chain behaves as a self avoiding walk with an effective step size per monomer 
given by 

<Mt = Y - (119) 

a (m + w) - 6 v ' 

From Table El it becomes clear that the effective monomer size changes from a value of 6 
for v = 0.00001 to a value of ~ 0.66 at the transition. The reason for the large value of 
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the effective monomer size at very small v is that the chains makes long jumps to take 
advantage of deep wells of the random potential, very much like anchored chains in a 
random potential that make sub-ballistic jumps [ 13j. This is also the reason why the 
number of monomers w in each well is somewhat less than L m , since a sufficient number 
of monomers need to be used for the connecting segments. For v 3> v c in the delocalized 
phase we expect the effective monomer size to be about 1, since the chain behaves almost 
like an ordinary self-avoiding chain with the random potential not playing any significant 
role. Thus the chain is expected to have its smallest size in the vicinity of the transition. 

It is important to notice that the discussion above is in the limit for very large L. If 
L < l/(2vG(V)), then in the localized phase when i;<» c the polymer will be confined 
to a single well and will appear compact, even though it will not remain so for large L. 
This may explain why in simulations that were done typically with L < 320 [ |2Ej the 
derealization transition appeared as a transition from a compact to a non-compact state 
of the polymer. 

We should also note that for an annealed random potential there is a transition from a 
collapsed state into an ordinary self-avoiding chain as v increases through the point v = g 
[ 111] ■ This is because the annealed free energy reads 

r™ L gL 2 vL 2 R 2 , , 

FiR ^W~-W + W + T- (120) 

For v < g the fourth term is negligible and the free energy is lowest when R — > (when 
L is large). For v > g the first term is negligible, and the radius of gyration grows like 

R ~ [v - g) l l 3 L 3/ \ (121) 

4.2. A self avoiding chain in a sea of hard obstacles 

We now turn to the case of hard obstacles. This case of was discussed at length in the 
previous section in the absence of a self-avoiding interaction. We will make use of the 
results applicable to three spatial dimensions. Three different behaviors were identified 
as a function of the system's volume. Region I is defined when the system's volume 
V < Vi ~ exp(ci/y / x), where c\ is a constant of order unity. Here < x < 1 is the average 
concentration of obstacles per site (total number of obstacles divided by total number of 
sites). We also assume that x is less than the percolation threshold (x c = 0.3116 for a 
cubic lattice), so sites occupied by obstacles don't percolate. In Region I we recall from 
the previous section that in the absence of a self-avoiding interaction, the free energy per 
site for a chain situated in a spherical region of volume R 3 in three dimensions is given 
by 

FilL = - \n(z) + 1/R 2 + x (122) 

where x is the actual concentration of obstacles in that region, whose minimal expected 
value in a system of total volume V is 



/ x In V 

x m ~x-y-^. (123) 

The binding energy per monomer inside the blob resulting from the lower concentration 
of obstacles in this region is given by x — x, since it is equal the entropy gain from a lower 
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concentration as compared to the average (background) concentration x. The chain is 
"sucked" towards regions with low concentration of obstacles since it can maximize its 
entropy there, and these regions of space act like the negative potential regions of the 
Gaussian random potential. Thus in order for the free energy per monomer to reflect 
correctly the binding energy of the chain inside the blob, both the constant x of the 
background and the constant term — ln(z), which is always there regardless of the chain's 
position, have to be subtracted. The relevant free energy per monomer situated in the 
blob (which is equal to minus the binding energy) is given by 



1 xlnV , s 

/' = ff-V-s- (124) 

This result coincides with Eq. ()23|) upon the substitution g — > x. Thus all the results of 
the previous section carry on to Region I with this simple substitution . 

Therefore we are going to discuss the situation when the system's volume is greater than 
V\ (Region II). In this case the many blob picture still holds, where the blobs are now 
situated in regions free of obstacles (with x — ) whose size is determined again by the 
distance Y of the jump which is also assumed to satisfy Y 3 > Vi ( an assumption which 
will be justified a posteriori) . In this case the farther the jump, it is more likely for the 
chain to find a larger space empty of obstacles, which will reduce further its confinement 
entropy w/R 2 and also the self-avoiding energy ( but there is a cost resulting from the 
connecting segments and the constraint of the total length being fixed). The blob size is 
given by R m ii = ^/G Q (the largest expected empty region in a volume Y 3 ) with [EI] 

Thus the free energy per monomer of a chain consisting of a number of blobs with con- 
necting parts, each of length Y, is given by 

f(w, m, Y) » - J— (Xl + ™+^- w{x - Go( ys )2) + vw 2 Go( Y*f) , (126) 
w + m \ m Y z Y A J 

where again the constant x, which is independent of m, w and Y, has been subtracted. 
This assures that the derealization transition again occurs at / = and not at / = x. 
The results of a numerical minimization of this free energy is displayed in Table El for the 
case of x = 0.1. 

The transition occurs between v = 0.041 and 0.042. Again we could find analytically 
almost the entire solution both for «<Ca; and at the transition. The solution is given in 
the Appendix. There we show that the transition occurs at v — 0.04142. We observe that 
as for the case of a Gaussian random potential the ratio w/m changes from ^> 1 to 1 
as v approaches the transition from below. The values of Y are seen to be consistent with 
the assumption Y 3 > Vi . We also checked that the free energy from Table II is lower 
than what one would obtain by constraining Y to be in Region I, i.e. Y 3 < V\. 

Finally for V > V2, where V2 — exp (x 2/(d+2) L d/(d+2) ) (Region III), we expect the be- 
havior of the chain to stay the same. This is because for jumps within a volume V2 the 
situation reverts to the previously discussed scenario, since the effective volume of interest 
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V 


Y 


m 


w 


f 


0.00001 


552 


2206 


72092 


-0.062 


0.0001 


127 


565 


9135 


-0.051 


0.001 


38 


206 


1241 


-0.034 


0.01 


23 


230 


207 


-0.0077 


0.02 


32 


536 


137 


-0.0019 


0.03 


59 


1737 


120 


-0.00035 


0.04 


194 


13915 


131 


-0.0000099 


0.041 


248 


21327 


134 


-0.0000023 


0.042 


357 


39310 


144 


0.0000024 



Table 3 

Free energy and parameters for a polymer in a sea of blockers when V > V%. 



which determines the statistics of the free spaces is of order Y 3 and we don't expect Y to 
be that large. 

An important point to note is that if one performs a simulation with strict self-avoiding- 
walks on a diluted lattice (with x < 1), one has v ~ 1, and hence one will always be in 
the delocalized phase and will not see any localization effects [l3~7~|. 

A few words are in order about the case of an annealed potential. This case has been 
already investigated in the literature [ I34j . and we will review it briefly. The free energy 
in the annealed case reads (for d = 3) 

F(R) * ± + xR* + ^. (127) 

The second term represents the entropy cost of a fluctuation in the density of obstacles 
that creates an appropriate spherical region of diameter ~ R. It was assumed that the 
chain occupies a spherical volume, or at least deviations from a spherical shape are not 
large [El]- In the case v = one obtains by minimizing the free energy that R ~ (L/x) 1 ^ 5 , 
a well known result. For v > the first term is irrelevant (and so is the "stretching term" 
of the form R 2 /L) and one finds that R ~ (f/x) 1 / 3 L 1 / 3 . In d dimension, the size scales like 
L x l d when v > 0, which is larger than the L 1/,( ^ +2 ) dependence in the v — case. There 
is no indication for a phase transition in these arguments, although some authors [ 13*4*] 
speculate that it breaks down for large v and a transition to a Flory L 3 /( d+2 ) dependence 
takes place. 

5. Conclusions 

In this chapter we have demonstrated the rich behavior of polymer chains embedded 
in a quenched random environment. As a starting point, we considered the problem of 
a Gaussian chain free to move in a random potential with short-ranged correlations. We 
derived the equilibrium conformation of the chain using a replica variational ansatz, and 
highlighted the crucial role of the system's volume. A mapping was established to that 
of a quantum particle in a random potential, and the phenomenon of localization was 
explained in terms of the dominance of localized tail states of the Schrodinger equation. 
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We also gave a physical interpretation of the 1-step replica-symmetry-breaking solution, 
and elucidated the connection with the statistics of localized tail states. Our concusions 
support the heuristic arrguments of Cates and Ball, but it starts with the microscopic 
model. 

We then proceeded to discuss the more realistic case of a chain embedded in a sea 
of hard obstacles. Here, we showed that the chain size exhibits a rich scaling behavior, 
which depends critically on the volume of the system. In particular, we showed that a 
medium of hard obstacles can be approximated as a Gaussian random potential only for 
small system sizes. For larger sizes a completely different scaling behavior emerges. 

Finally we considered the case of a polymer with self-avoiding (excluded volume) in- 
teractions. In this case it was found that when disorder is present, the polymer attains 
the shape like that of a pearl necklace, with blobs connected by straight segments. Us- 
ing Flory type free energy arguments we analyzed the statistics of these conformational 
shapes, and showed the existence of localization-delocalization transition as a function of 
the strength of the self-avoiding interaction. 

The work described in this chapter is concerned with static (equilibrium) properties 
of polymers in random media. There is a lot of theoretical work still to be done related 
to the dynamics, and especially nonequilibrium properties of polymers in random media. 
This is also of practical importance, for example for the separation of chain of different 
length or mass like DNA molecules under the effect of an applied force when embedded 
in a random medium like a gel [ 38 . 
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